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0 . ABSTRACT  | J\  \ | 

Models  are  given  for  sequences  of  correlated  exponential 
interarrival  and  service  times  for  a single  server  queue.  These 
multivariate  exponential  models  are  formed  as  probabilistic  linear 
combinations  of  sequences  of  independent  exponential  random  vari- 
ables and  are  easy  to  generate  on  a computer.  Limiting  results 
for  customer  waiting  time  under  heavy  traffic  conditions  are 
obtained  for  these  queues.  Heavy  traffic  results  are  useful  for 
analyzing  the  effect  of  correlated  interarrival  and  service  times 
in  queues  on  such  quantities  as  queue  length  and  customer  waiting 
time.  They  can  also  be  used  to  check  simulation  results. 


Support  from  NSF  grarvt»-^NG- 77^-090 20  and  ENG-79-01438  and 
ONR  contract  NR-42-284  isbae^tefully  acknowledged. 


1 . INTRODUCTION 


Much  of  the  work  in  queueing  theory  deals  with  models  in 

which  assumptions  of  independent  service  and  interarrival  times 

are  made.  One  way  to  examine  the  effect  of  correlated  interarrival 

and  service  times  on  measures  of  performance  in  queues  is  to 

consider  models  in  which  the  interarrival  and  service  times  are 

correlated  stationary  sequences  of  random  variables  having  given 

marginal  distributions  and  compare  results  for  these  models  with 

those  for  which  the  interarrival  and  service  times  are  independent 

and  identically  distributed  with  the  same  marginal  distributions. 

In  [3]  Jacobs  and  Lewis  introduced  a scheme  for  generating 

sequences  of  dependent  exponential  random  variables.  The 

EARMA  (1,1)  (exponential  mixed  autoregressive  moving  average  with 

both  autoregression  and  moving  average  of  order  1)  model  is  defined 

as  follows.  Let  f E } be  a sequence  of  independent  exponential 

random  variables  with  positive  finite  mean  v-1.  Let  {J  } and 

n 

{ Kn } be  independent  sequences  of  independent  (0,l}-random  variables 
with  P{Jn=l)  = 1-6  and  P{Kn=l}  = 1-p  where  0 <_  6 £ 1 
and  0 ^ p < 1 are  given  constants.  We  define 

<1*1)  Xn  = BEn  + JnAn-l  ' n = 1 

where 

U-2)  An  = pAn_1  + KnEn  , n = 1,2,...  . 

We  assume  that  Ag  is  an  independent  exponential  random  variable 
with  mean  v Under  this  assumption  (Xn)  is  a stationary 

sequence  of  dependent  exponential  random  variables  with  mean  \>  E . 
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The  dependence  is  not  Markovian  in  general  and  as  was  shown  in  [3] 


(1.3)  Corr(XrX1+n)  = pn_1  (1-6  ) t 6(l-p ) + (l-3)p  ] . 

Note  that  if  8=1  or  8=0  and  p = 0,  then  (Xn)  i6  a 
sequence  of  independent  exponential  random  variables. 

The  EARMA (1,1)  process  is  formed  as  a probabilistic  linear 
combination  of  independent  random  variables.  This  technique  can 
be  extended  to  form  sequences  of  multivariate  exponential  random 
variables.  This  idea  was  briefly  explored  in  Lewis  and  Shedler  [7]. 
These  models  have  the  advantage  that  the  marginal  distribution  and 
correlation  structure  of  the  sequence  are  specified  separately. 

In  section  2 we  will  give  several  ways  in  which  to  generate 
sequences  of  dependent  multivariate  exponential  random  variables 
from  probabilistic  linear  combinations  of  independent  exponential® . 
These  stationary  multivariate  EARMA  sequences  of  dependent  exponentials 
can  be  used  as  interarrival  and  service  times  in  queues.  Results 
for  measures  of  performance  for  these  queues  can  then  be  compared 
to  queues  with  stationary  independent  sequences  of  interarriva’l 
and  service  times  having  the  same  marginal  (exponential)  distri- 
butions. Since  the  correlation  structure  in  the  multivariate 
EARMA  sequences  is  specified  by  parameters  independent  of  the 
marginal  distribution,  any  difference  between  results  for  queues 
with  these  multivariate  interarrival  and  service  sequences  and 
the  corresponding  queues  with  independent  interarrival  and  service 
sequences  having  the  same  marginal  distributions  will  be  due  to 
the  correlation  structure. 


Suppose  that  the  interarrival  and  service  times  of  a single 
server  queue  form  stationary  multivariate  EARMA  sequences 


'It 

of  dependent  exponential  random  variables.  Results  of  Loynes  [9] 
give  conditions  for  the  existence  of  limiting  distributions  for 
the  waiting  time  of  a customer  and  the  queue  length.  However,  it 
seems  to  be  difficult  to  obtain  exact  results  even  in  the  simplest 
queues  with  EARMA  interarrival  and  service  times  (cf.  Jacobs  [4]). 

As  a result  we  turn  our  attention  to  the  possibilities 
of  obtaining  limiting  results  of  the  heavy  traffic  type  for 
queues  with  stationary  multivariate  EARMA  sequences  of  dependent 
interarrival  and  service  times.  Results  of  this  kind  will  indi- 
cate how  dependence  in  queues  introduced  by  the  EARMA  scheme 
affects  quantities  such  as  customer  waiting  time  and  queue  length, 
at  least  under  heavy  traffic  conditions.  These  limiting  results 
can  also  be  used  to  check  simulations. 

In  this  paper  we  will  show  that  Kingman's  [5]  result 
for  the  equilibrium  waiting  time  W in  heavy  traffic  holds  for 
queues  having  multivariate  EARMA  interarrival  and  service  times; 
that  is,  we  will  show  under  conditions  detailed  in  Section  3,  that 
if  the  traffic  intensity  Ay  ^ of  the  queue  is  close  to  1,  then 
the  distribution  of  (1  - Ay  is  approximately  exponential. 

The  mean  of  this  exponential  distribution  will  in  general  be 
different  from  the  case  of  independent  interarrival  and  service 
times.  Positive  cross-correlation  between  the  interarrival  and 
service  times  will  tend  to  decrease  the  mean;  while  positive 
correlation  within  the  interarrival  and  service  time  sequences 
will  tend  to  increase  it. 
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In  Section  3 we  will  state  the  heavy  traffic  result  more 
precisely  and  present  several  examples  of  its  application-  Its 
formulation  for  multivariate  EARMA  queues  is  very  dependent  on 
the  fact  that  the  marginal  distribution  and  correlation  structure 
are  specified  independently  of  one  another  in  EARMA  formulations- 
In  Section  4 we  will  give  the  proof  of  the  heavy  traffic  result. 


(2.1)  EXAMPLE.  Let  (Xn>  be  an  EARMA (1,1)  process  with 

positive  finite  mean  and  parameters  0 < 6^  < 1 and 

0 £ Pi  < 1 ' ^et  (Sn>  1:36  an  process  independent  of 

{Xn>  with  positive  finite  mean  y-1  and  parameters  0 £ 32  £ 1 
and  0 £ p2  <1.  Note  that  if  $^  = 1,  $2  = 1,  or  3^=1, 

$2  = 0 , p2  = 0,  or  3^  = 0,  Pjl  = Of  32  = ^ • or  6-^  = 0 , = 0 

$2  ~ 0 f P2  = 0 , then  the  queue  reduces  to  an  M/M/1  queue . 

Other  special  cases  include  queues  with  Poisson  arrivals  and 
EARMA (1,1)  service  times  independent  of  the  arrival  process, 
and  queues  with  EARMA (1,1)  interarrival  times  and  independent 
exponential  service  times. 

The  sequence  {(Xn,Sn)}  is  stationary  and  {Xn>  and 
{Sn}  are  independent.  However  the  xn's  (respectively  Sn's) 
are  positively  correlated  with 


(2.2)  Cov(Xn,Xn+k)  = X_2pJ"1(l-ei) (61(l-p1)  + (1-61)P1] 
(respectively 

(2.3)  Cov  (sn'Sn+K  ' = (l-e2>  IB2  (1-P2)  + d-B2>  P2J ) 

from  (1.2). 

1 2 

Let  Fm  (respectively  F^)  be  the  o-algebra  generated 

by  X1,...,Xm  (respectively  S1,...,Sm)  and  G*+k  (respectively 
2 

Gm+k)' be  the  o-algebra  generated  by  Xm+k,  xm+)c+x***-  (respectively 

Sm+k*  Sm+k+l'  •••>•  L2(Em>  'respectively  L2(S*+k>> 

denote  the  collection  of  real-valued  functions  having  finite 
second  moment  that  are  measurable  with  respect  to  (respective! 

Sj,+k>-  1 ■ 1'2‘  I£  fi  4 *nd  % * l.2tgi+k)-  2 * 2’2' 

then  there  is  a kg  independent  of  f^,  g^,  i = 1,  2,  such  that 
for  k 2 kQ 

(2.4)  |E[  (f1f2)  (gxg2)  ] - E[f1f2J  Kigali 

= |E[figil  E[f2g2]  - E[f1J  E[f23  E[gx]  E[g2]| 

» |E[figi]  {E[f2g2]  - E [f 2 ] E(g2]} 

+ E [f 2 ] Etg2]  {E[figil  - E[fx]  Efg^Jl 

6 

— ™ -r  ;x7  ;r; 

, y„  -■  tea- 
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(2.5) 


< E [ | fl9l  I ] 5pk/2  E [ f 2 ] 1/2  E[g2]1/2 

+ E [ | f 2 | ) E [ | g 2 I 1 5pk/2  E[f2]1/2  E[g2]1/2 

< 5[pk/2  + Pk/2]  e-[  (f1f2>  2]1/2  E [ (gig2)  2] 1/2 


where  (2.4)  follows  from  the  independent  of  and 

and  (2.5)  follows  from  (6.5)  of  Jacobs  and  Lewis  [3].  Hence, 

if  f (respectively  g)  has  finite  second  moment  and  is  measur- 

1 2 

able  with  respect  to  the  product  o-algebra  £m  x Em  (respectively 
-m+k  X -m+k* ' then  for  k sufficiently  large 

(2.6)  |E[fg]  - E [f ] E[g]|  < 5E[f2]1/2  E[g2]1/2  [p\/2  + 

Hence,  the  sequence  {(X  ,S  )}  is  asymptotically  uncorrelated 
^ n n 

in  the  sense  of  Rosenblatt  [10]  and  is  ^-mixing  in  the  sense  of 
3illingsley  [1]. 

(2.7)  EXAMPLE.  In  Lewis  and  Shedler  [7]  one  of  the  schemes 

that  is  proposed  for  cross-correlating  the  arrival  process  and 

the  service  times  is  the  following.  Assume  {Xn}  and  (En} 

are  independent  sequences  of  independent  exponential  random 

variables  with  positive  finite  means  A ^ and  y * respectively. 

The  nth  interarrival  time  is  assumed  to  be  X and  hence  the 

n 


arrival  process  is  Poisson  with  rate  A.  Let  the  nth  service  time 


E 


where 


(2.9) 


Bn  = 


PB 


n-1 


+ K X 
n n 


and  8,  p,  (J„}  and  (K  } are  as  in  (1.1)  and  (1.2)  so  that 
n n 

(Sn>  and  {Xn>  are  now  cross-correlated  sequences.  We  will 
assume  that  BQ  has  an  exponential  distribution  with  mean  A 

Since  En  and  Bn  are  independent,  { Sn } is  a stationary 
sequence  of  exponential  random  variables  with  mean  y~* . If 
8=1,  then  the  queue  reduces  to  an  M/M/1  queue.  The  Sn's 
can  be  shown  to  be  positively  correlated  with 

(2.10)  Cov(Sn,  Sn+k)  = y-2 (1-8) 2 Pk,  k = 1,2,...  . 

Further,  the  nth  service  time  is  positively  correlated  with 

the  interarrival  times  X . X^  , ...  , X.  with 

n n-i  l 


(2.11) 


Cov  (S 


n' 


xn_k)  = (Ay)_1(l-8)  (l-p)pk 


for  k = 0 , . . . ,n-l . 

Note  that  (X  ,Sn,...,  Xn+k,Sn+Jc)  is  a function  of 

(Bn-l' VEn'Jn' VXn+l'En+l'Jn+l'Kn+l'  * ** ' Xn+k,En+k,Jn+k'Kntk) 
Hence,  since  the  sequence  is  stationary  so  is  the 

sequence  { ( Xn , Sn ) ) - Further,  by  arguments  similar  to  those 

used  in  obtaining  (6.5)  of  Jacobs  and  Lewis  [3],  if  f and 

g are  as  in  (2.6),  then  for  k sufficiently  large 


(2.12) 


I E [fg]  - E [f ] E[g]|  < 5pk/2  E [f2] 1/2  E[g2]1/2 


so  that  the  sequence  {(Xn,Sn)}  is  ^-mixing. 
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There  are  many  ways  to  cross-correlate  the  interarrival 
and  service  times  using  the  EARMA  scheme  beyond  those  given 
above  and  in  Lewis  and  Shedler  [7].  We  will  present  one  more, 
somewhat  more  general  than  that  of  Example  (2.7)  in  which  the 
service  times  were  coupled  back  over  the  previous  interarrival 
times . 


(2.13)  EXAMPLE.  Let  {En> , {En>,  and  {Cn>  be  independent 
sequences  of  independent  exponential  random  variables  having 
respective  positive  finite  means  A-1,  y*"^ , and  1.  Let 


(2.14) 


X = BE  + J U_1A  ) 
n n n n 


(2.15) 


Sn  * + 


with 


(2.16) 


A = pA  , + K C 
n H n-1  n n 


where  6,  p,  {Jn>  and  {Kn}  are  as  in  (1.1)  and  (1.2)  and  (In> 
is  a sequence  of  independent  binary  random  variables  independent 
of  everything  with  P{In=l}  = 1 - P{In=0}  = (1-y)  for  some 


fixed  constant  0 £ y £ 1.  Again  we  assume  that  «.q 


Art  is  an. 

independent  random  variable  having  an  exponential  distribution 


with  mean  1 . Then  {(E,E,A,I,J)}  is  a stationary 

n n n n n 

sequence  of  random  variables  and,  hence,  the  same  is  true  of 

(<VSn)K 


' MB  m'to  trtii  <MI  m 'I  


«■'■  ■■  ■ ■ . 


lllUflftiiififlt'T  |- 


If  Y = 1.  8=1  or  Y = 1/  8=0,  p = 0 or  y = 0, 

8=1,  p = 0,  then  the  queue  reduces  to  an  M/M/1  queue.  If 
8=1,  p = 0 then  the  queue  reduces  to  a special  case  of  Example 
(2.7).  Straightforward  manipulation  shows  that  in  general  the 
interarrival  times  Xn  are  positively  correlated  with 

(2.17)  Cov(Xn,  Xn+k)  = A~2(l-8)2pk  , k = 1,2, 

The  service  times  Sn  are  also  positively  correlated  with 

(2.18)  Cov(Sn,  Sn+k)  = y~2(l-Y)2pk,  k = 1,2,...  . 

The  interarrival  times  and  service  times  are  positively  cross- 
correlated  with 

(2.19)  Cov(Xn,  Sn+k)  = (1-6)  (1-Y)  (Xy)-1pk,  k = 0,1,...  ; 
and 

1 ' 

■ i 

(2.20)  Cov(Xn+k,  Sn)  = (1-6)  (1-y)  avi)_1pk,  k = 0,1,...  . 

Finally,  by  the  arguments  used  in  obtaining  (6.5)  in 
Jacobs  and  Lewis  [3],  if  f and  g are  as  in  (2.6),  then 

i 

(for  k sufficiently  large 

(2.21)  | E [ fg]  - E [f ] E [g] I < 5pk/2  E [f 2] 1/2  E[g2]1/2  . 

so  that  the  sequence  { ( Xn , sn ) } is  (((-mixing. 
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3.  A HEAVY  TRAFFIC  RESULT  FOR  EARMA  QUEUES 

Let  {Xn>  (respectively  (Sn>)  be  the  sequence  of 

interarrival  times  (respectively  service  times)  for  a single 

server  queue  with  first- in-first-out  service  discipline. 

Put  U = S„  - X and  assume  that  {U  } is  a stationary 
n n n n 

sequence  with  E[Unl  < 0*  Let  Wn  be  the  waiting  time  of  the 

nth  customer  and  W = supn  > ^ l£=l  Uk  (with  the  convention 

that  = 0).  Loynes  [9]  has  shown  that  if  the  sequence 

{U  } is  metrically  transitive,  then  the  distribution  of  W 
n n 

tends  in  the  limit  to  the  distribution  function  of  W.  Since 
it  is  usually  difficult  to  obtain  the  distribution  of  W 
analytically,  particularly  when  the  service  and  interarrival 
times  are  dependent,  we  will  obtain  a heavy  traffic  limiting 
result  for  it.  In  this  section  we  will  give  the  statement 
of  the  result  as  well  as  some  examples  of  its  use  leaving  the 

proof  to  the  next  section. 

For  each  a > 0 let  Qa  be  a single  server  queue 

with  FIFO  queueing  discipline  having  stationary  sequences  of 
interarrival  times  (Xn(a)}  with  positive  finite  mean  A (a)-'1' 
and  service  times  (Sn(a)}  with  positive  finite  mean  y(a)-^ 
such  that,  if  Un (a)  = Sn  (a)  - Xn(a),  then  (Un(a)}  is  a 
stationary  sequence  of  random  variables.  Let 

(3.1)  E[U  (a)]  = v(a)  < 0;  Var[U  (a)]  = s2  (a)  ; 
n n 
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mut 


v.  . 


(3.3)  d (a)  = | v (a)  | /a  (a)  and  n(a)  = d(a)~^  . 

We  will  assume  the  following: 

(3.4)  ASSUMPTION.  There  exists  a distribution  function  G so 

that  lim  P{[U  (a)  - v(a)]2  < x}  = G(x)  for  all  continuity 

• 2 2 
points  x of  G;  lima  ^ ^ v (a)  = 0 ; and  0 < lima+o(s  (a)  = s < ». 

(3.5)  ASSUMPTION.  o < lim  a2( a)  = a2  < ~ . 

a -*■  00 

(3.6)  ASSUMPTION.  If  f (respectively  g)  is  a function  with 
finite  second  moment  measurable  with  respect  to  the  a-algebra 
generated  by  ^ (a)  , . . . ,Um(a)  (respectively  Um+k(a),  um+k+1  <<*>  » • • • ) » 
then 

(3.7)  |E[fg]  - E[f]  Elg]|  < <)>(k)  E[f2)1/2  E[g2]1/2 

for  some  decreasing  function  <J>  bounded  above  by  1 not  depending 
on  a satisfying 


(3.8)  l <p  (n)  < ~ . 

n=l 

Assumption  (3.4)  concerns  the  convergence  of  the  distribution  of  U2 (a) 

2 n 
— (o)  (ot)  J and  its  first  two  moments  as  a Assumption 

(3.6)  is  roughly  a condition  that  the  sequences  (u  (a) , n = l,2,...} 
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be  4>-mixing  uniformly  over  all  a.  This  assumption  guarantees 
that  the  series  in  (3.2)  converges  and 

lim  - E[{  Y (U.  (a)  - v (a) ) ) 2 ] = a2(a)  . 
n ~ n k=l 

Assumption  (3.5)  concerns  the  convergence  of  this  normalized 
variance  term  as  a 

Put 

, [n(a)t] 

Z (t)  = d (a)  o (a) ~ I Uk(a),  t >_  0 
a k=l 

where  [x]  is  the  largest  integer  less  than  or  equal  to  x. 

Let  3D([0,«))  be  the  collection  of  real-valued,  right  continuous 
functions  on  [0,«>)  which  have  left  limits  everywhere  endowed 
with  the  metric  given  in  Lindvall  [8]. 


(3.9)  THEOREM.  If  (3.4) -(3.6)  are  satisfied,  then  the  process 
Za  = (Za(t);t  > 0}  converges  weakly  in  D([0,°°))  to  a 
Brownian  motion  with  negative  unit  drift. 

The  proof  of  this  result  will  be  given  in  the  next  section. 
Let 

n 

W(a)  = sup  l U.  (a)  . 
n > 0 k=l  K 

Then  it  follows  from  Theorem  (3.9)  that 
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lim  w(a)  < x|»  = lim  P(  sup  Z (t)  < x} 

a oo  lo(a)  J a -*■  00  t^O 

= 1 - e xp  { - 2 x } 

where  the  last  equality  follows  from  the  continuous  mapping 
theorem  (cf.  Billingsley  [1]).  Thus  we  have  the  following 

result. 

(3.10)  THEOREM.  If  (3.4)— (3.6)  hold,  then 

^ 2 

lim  p{ | v (a) | W(a)  < xJ  - 1 - exp(-2a_  x>  • 

a 00 

Hence,  if  the  traffic  intensity  Ay  1 of  a queue  is  less 
than  but  close  to  1 and  (Un>  is  a stationary  sequence  satisfy- 
ing (3.7)  and  (3.8),  then  the  distribution  of  W is  approximately 

1 2 -1 

exponential  with  mean  - j a ElU^l  where 

_ 00 
o = Var  [U1  ] + 2 \ CovOJ^U^)  . 

n=l 

Thus  positive  cross  correlation  between  the  interarrival  and  service 
times  will  tend  to  decrease  the  average  mean  waiting  time;  while 
positive  correlation  within  the  interarrival  time  or  service 
time  sequence  will  tend  to  increase  the  mean  waiting  time. 
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(2.11)  EXAMPLE.  For  each  a > 0,  let  Qq  be  a queue  as  in 
Example  (2.1)  with  {Xn(ot)}  (respectively  {Sn(a)})  forming 
an  EARMA (1,1)  process  with  parameters  X(a),  B^  and 
(respectively  y(a),  B2  and  p2) . Assume  X(a)  < y(a)  for 
each  a,  but  lim  X(a)  = lim  ^ y (ot)  = m;  then  hypotheses 
(3.4) -(3.6)  are  satisfied  and  Theorem  (3.10)  holds  with 

a2  = 2m~2 (1  + (l-P]^)-1  (l-B^  [(l-p1)B1  + (l-B^p^ 

+ (1-P2)-1  (1-B2)  [(1-p2)B2  + d-62)p2]>  . 

Hence,  if  the  traffic  intensity  Xy-*  of  a queue  with  independent 

EARMA (1,1)  services  and  interarrival  times  is  less  than  but 

close  to  1,  then  the  distribution  of  the  equilibrium  waiting 

1 2 

time  W is  approximately  exponential  with  mean  - j a E[U^] 

-1  -1  2 

where  E[U^]  = y - X and  a is  given  by 

a2  = X"2  {1  + 2(l-pir1  d-Bx)  [(l-p1)B1  + (1-B1)P2]} 

+ y-2{  1 + 2(1— p2 ) ~ 1 (1-B2)  [(l-p2)32  + (1-B2)P2)}  . 

Thus  the  dependence  in  this  queue  tends  to  increase  the  mean 
customer  waiting  time  at  least  in  close  to  heavy  traffic 
conditions . 
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(3.12)  EXAMPLE.  For  each  a > 0,  let  be  a queue  as  in 

Example  (2.7)  with  parameters  E[Xn(a)J  = A (a)-1,  3,  and  p 
and  E[Sn(a)]  = y(a)"1.  Again  we  will  assume  A (a)  < p(a) 
for  each  a but  lima+M  A(n)  = A = lima  ^ ^ p(a)  = p = m.  Hence 
(3. 4) - (3 .6)  are  satisfied  and  Theorem  (3.10)  holds  with 


°2  - A-2  + p"2  - 2(1-3)  (1  - p)  (Ap)_1 

+ 2p(l-p)_1  [ (1-3)  2p-2  - (1-3)  (1  - p)  (Ap)"1] 


-2  . -2 


= A + p [1  + 2 (1-3)  p(l-p)’1]  - 2(Ap)“i  (1-3) 


= 2m-2  [1  + (1-3) 2 p(l-p)'1  - (1-3)]  . 


2 2 

The  quantity  a - aind  can  be  either  positive  or  negative  where 
°2nd  = A 2 + p 2 = 2m  2.  It  is  negative  if  p < (2-3) 
and  nonnegative  otherwise.  Hence,  if  p is  sufficiently  less 
than  (2-3)  1,  then  the  average  customer  waiting  time  in  the 
queue  will  tend  to  be  less  than  the  waiting  time  in  the  inde- 
pendent case  at  least  in  close  to  heavy  traffic  conditions. 

It  will  tend  to  be  greater  if  p is  sufficiently  greater  than 
(2-3) -1 . 


(3.13)  EXAMPLE.  For  each  a > 0,  let  be  a queue  as  in 

Example  (2.13)  with  parameters  A (a),  p(a),  p,  3, and  y- 


Assume  lim  ^ A (a)  = A = lim 

a-*-  °°  or* 


p = m as  before. 


Then  (3. 4) -(3. 6)  are  satisfied  and  Theorem  (3.10)  holds  with 
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o2  = p"2  + X-2  - 2 (Xp) -1  (1-3) (1-Y) 

+ 2p(l-p)"1{y-1  (1-y)  -X-1 (1-3) )2 

= 2m~2  [l-(l-3)  (1-y)  +p(1-p)-1{(S-y)}2]  . 


2 2 

The  difference  a - can  be  either  positive  or  negative; 

2 -2  -2  -2 

(o^n<j  * X + p = 2m  ) . The  difference  will  be  negative  if 

2 

p < (1-8) (1-y) /[  (1-3)  (1-y)  + (B-y)  1 and  nonnegative  otherwise. 

Hence  under  close  to  heavy  traffic  conditions,  the  average 

customer  waiting  time  will  tend  to  be  less  than  the  independent 

2 

case  if  (1-8)  (1~y) / 1 (1-6) (1-y)  + (8-y)  1 is  sufficiently  greater 
than  p.  It  will  be  greater  than  the  independent  case  if  the 
same  quantity  is  sufficiently  less  than  p. 
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4.  PROOF  OF  THE  MAIN  RESULT 

In  this  section  we  will  give  a proof  of  Theorem  (3.9). 
We  will  assume  that  (3.4) -(3.6)  are  satisfied.  Let 


y (t)  - did 

Yalt)  - 3TST 


[Uk(ot)  - v (a)  ] 


By  Theorem  3'  of  Lindvall  [8]  and  the  continuous  mapping  theorem. 
Theorem  (3.9)  will  be  proved  if  for  each  a > 0 we  can  show 
the  following  result. 

(4.1)  THEOREM.  The  process  {Ya(t);  0 £ t £ a}  converges 
weakly  on  3D(fO,a])  to  Brownian  motion. 


0D([O,a])  is  the  space  of  real -valued  right  continuous  functions 
on  [0,a]  with  left-hand  limits  endowed  with  the  Skorohod  ^ 
topology . ) 


We  will  prove  Theorem  (4.1)  for  a = 1 by  a series  of 
lemmas.  The  proof  for  arbitrary  a > 0 is  similar. 

(4.2)  LEMMA.  The  process  ( t) ; 0 £ t < 1}  has  asymptotically 
independent  increments. 

Proof.  Let  0 £ s1  £ t1  < s2  £ t2  < •••  < sr  £ t £ 1.  Let 
ci  be  the  event  (Y^tj.)  - YQ(si)  £ xi>  for  real  X;L . Then 
is  an  event  in  the  o-algebra  generated  by 
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If  6 is  the  smallest  difference 


U[n(a)s..]+1'  ' U[n(a)t.)‘ 

Sj^  - then  [n(a)si]  + 1 - [n(a)ti_1]  >.  [n(a)6].  By  (3.7) 

I P ( 0 C.  ) - n P(C.)  I < r*([n(a)«])  . 
i=l  1 i=l  1 

Since  6 > 0,  r$([n(a)6])  tends  to  zero  as  a -*•  « and  the 
result  follows. 

In  the  following  results,  let  (a)  = IL(a)  - v (a) 

and  V (a)  = Tf  , D.  (a) . By  (3.4)  — (3.6)  the  following  result 
m 1=1  l 

holds . 


(4.3)  LEMMA. 


a)  £E[Vk(c02 *]  <2  ^ |E[D1(a)  D1+j(a)]|; 

b)  lim  n (a) -1  E[V  . .(a)2]  = a2  ; 

a > ' ' 


c) 


If  is  bounded  by 

where  K,  depends  on 

<p 


C,  then  Etvn(a)  (a)4]  <_  K^C4  n(ct)2 
<p  alone . 


(4.4)  LEMMA.  The  sequence  of  random  variables 
-1  2 

(n(a)  (vn(a)  (“)  ) » a >.  0)  is  uniformly  integrable. 

2 

Proof.  By  (3.4),  {D^(a);oO  0}  is  uniformly  integrable.  The 

remainder  of  the  proof  is  very  similar  to  that  of  (20.48)  on 
page  176  of  Billingsley  [1]  and  will  be  omitted. 


? 


i 
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(4.5)  LEMMA 


such  that  n(a)  > n, 


P{  max  |v.  (a)  | > y/n  (a)  } < 
i < n (a)  1 


such  that  n(a)  > m 


then 


By  the  previous  lemma  we  can  choose 


P{  max  |V.  (a)  | >_  3y  /n  (a)  } 
i <_n(a) 

< P(  |Vn(a)  (a)  | > y OT7-) 

+ [n(a)  +p(a)]  p{Cp(a)(a)  1 Y 1/11  (a)  } 

n (a) -p (a) -1  . , 

+ .I  P(Ai(a)  [P{  IVn(Qt)  (a)  - vi+p(a)  <a>l>Y  /n  (a)  }+<Mp  (a)  ] 

< 2eY_2+<t»(p  (a)  ) + [n(a)  +p(a)]  p{C*(a)(a)  y •/nTaT> 

by  (4.7)  and  the  fact  that  the  A^(a)  are  disjoint.  Hence  by 
(4.6)  and  (3.8)  for  a large 

P{  max  |v.  (a)  | 3y  /n  (a)  ) <_  3ey  2 . 
i <_ n (a)  1 

PROOF  OF  THEOREM  (4.1).  For  a = 1:  By  Lemma  (4.5)  the 

sequence  {Y  } is  tight  on  3D [0,1],  Further  lira  ^ E [Y  (t)]  =0 

and 


by  (3.5).  Since  Ya(t)2  = [a2(a)  n(a)]-1  V[n(a)t](a)2  , 

2 

{ Y^ ( t ) } is  uniformly  integrable  for  each  t by  Lemma  (4.4). 

By  Lemma  (4.2)  {Ya(t);  0 <_  t <_  1}  has  asymptotically  independent 

increments.  The  result  now  follows  from  Theorem  (19.2)  of 
Billingsley  [1]. 


lim 


a 


E[Ya(t) 


= lim 
a -*■ 


(/n  (a)  a (a)  ) 


-2 


[n  (a)  t] 

l D 

k=l  k 


(a) 


(4.7)  p{  I Vk  (a)  (“>  I > Y ^ ^ < eY~^ 

for  all  a.  Let 

A.  (a)  = {max  |v.  (a)  | < 3y  /n  (a)  < | V . (ot)  | > ; 
j < i J 1 

then 

P{  max  |v. (a) | > 3y  /n (a) } 
i £ n ( a ) 

£ p{  lvn(a)  (a)  | >_  y /n  (a)  } 
n (a) -1 

+ J P(Ai(cx)  n (IVn(a)  -Vi(a)|  > 2y  ) . 

The  summation  is  bounded  above  by 


n (a)  -p  (a)  -1 


P{  |vi(a)  - vi+p(a)  («)  I > y (a)  } 


n (a)  -p  (a)  -1 

+ P(A.(a)  n {IVn(a)(a)  ~ V.+p(a) (a) | > y /HJaT} ) 

n (a)  -1  

+ l p^Vnfntl*a)  - v.  (a)  | > y ^n  (a)  } 

i=n  (a)  -p(a)  n(oU  1 


Each  term  in  the  first  and  third  of  these  sums  is  bounded  above 
by  P{Cp(a)  (a)  >_  y /n  (a) } . Since  A^(a)  is  measurable  with 
respect  to  the  o-algebra  generated  by  U^ot),...,  U^a)  by  (3.7) 
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5 . EXTENSIONS 


This  paper  presented  a heavy  traffic  result  for  the 
limiting  distribution  of  the  customer  waiting  time  in  single 
server  queues  with  dependent  interarrival  and  service  times. 
There  are  many  different  ways  in  which  to  model  dependent 
interarrival  and  service  times  in  queues  using  the  EARMA  scheme . 
Some  other  ways  are  given  in  Lewis  and  Shedler  [7],  Properties 
of  these  multivariate  exponential  models  will  be  given 
elsewhere . 

Work  is  currently  going  on  in  the  simulation  of  these 
queues.  Some  simulation  results  have  been  reported  for  the 
queue  of  Example  (2.7)  with  6 = 0.5,  y = 4 and  various 
values  of  P and  A in  Boonsong  [2].  They  seem  to  indicate 
that  the  heavy  traffic  result  approximates  the  average 
customer  waiting  time  fairly  well  for  a traffic  intensity  of 
.9  and  quite  well  for  a traffic  intensity  of  .95.  However, 
the  distribution  of  W in  a simulation  with  p = 0.25  and 
traffic  intensity  .995  shows  a large  amount  of  underdispersion 
relative  to  an  exponential  distribution.  Thus,  it  seems  that 
the  convergence  of  the  waiting  time  distribution  is  very  slow. 

It  is  expected  that  further  simulations  will  give 
more  indication  of  the  validity  of  using  the  heavy  traffic 
result  to  approximate  the  average  customer  waiting  time,  or 
quantiles  of  the  waiting  time  distribution  under  less  than 
heavy  traffic  conditions.  These  results  will  be  reported 
on  elsewhere. 
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